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1. Introduction 

Quantum Chromodynamics (QCD) — the underlying theory of the strong interactions — when 
combined with the electroweak interactions, is responsible for all of nuclear physics, from the 
structure of light nuclei to the most complex nuclear reactions. Presently lattice QCD is the only 
known method which allows one to compute properties of hadrons and nuclei from first principles. 
While lattice QCD provides highly-accurate predictions and postdictions of hadronic properties 
Uke masses and decay constants, the study of hadronic interactions, like two-body scattering and 
three-body interactions is still in a development phase. The NPLQCD collaboration has calculated 
the interactions between many of the lowest-lying baryons and mesons, including 7r+7r+, 7c'+7r+, 
K^K^ , meson-baryon, nucleon-nucleon, and hyperon-nucleon scattering as well as multi-meson 
systems (much of this work is reviewed in Ref. It is noteworthy that the 71+71+ scattering 
length has been calculated to percent-level accuracy at the physical point [^. Part of the reason 
behind this accuracy is that QCD correlation functions involving Goldstone bosons are subject to 
powerful chiral symmetry constraints. Since current lattice calculations are carried out at unphys- 
ical quark masses, these constraints play an essential role in extrapolating the lattice data to the 
physical quark masses, as well as to the infinite volume, and continuum limits. Chiral perturbation 
theory {x-?T) is the optimal method for implementing QCD constraints due to chiral symmetry, 
and in essence, provides an expansion of low-energy S -matrix elements in quark masses and powers 
of momentum 

In contrast to the purely mesonic sector, recent studies of baryon-baryon interactions — the 
paradigmatic nuclear physics process — have demonstrated the fundamental difficulty faced in 
making predictions for baryons and their interactions [Q, |5|]. Unlike the case with mesons, correla- 
tion functions involving baryons suffer an exponential degradation of signal/noise at large times |^]. 
Furthermore, while baryon interactions are constrained by QCD symmetries like chiral symmetry, 
the constraints are not nearly as powerful as when there is at least one pion or kaon in the initial or 
final state. For instance, there is no expectation that the baryon-baryon scattering lengths vanish in 
the chiral limit as they do in the purely mesonic sector. To make matters worse, in s-wave nucleon- 
nucleon scattering, the interactions are enhanced due to the close proximity of a non-trivial fixed 
point of the renormalization group, which drives the scattering lengths to infinity, thus rendering 
the effective field theory description of the interaction non-perturbative |^. 

In this proceedings I will review the most recent progress made in the description of hadronic 
scattering and interactions. Given the contrast in difficulty between the purely mesonic and purely 
baryonic sectors described above, it is clearly of great interest to perform a lattice QCD investiga- 
tion of the simplest scattering process involving at least one baryon: meson-baryon scattering. The 
first fully-dynamical study of meson-baryon scattering has been carried out recently in Ref. [||] and 
will be reviewed. I will also review other recent work [§, [T^ by the NPLQCD collaboration, which 
has performed very-high statistics calculations of few-body baryon systems that offer new insights 
into the signal/noise problem, and allow a first glimpse of the baryon number B = 3 sector. 

2. Signal to Noise Estimates 

As is well known [^], very general field-theoretic arguments allow a robust estimate of the noise to 
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signal ratio of hadronic correlation functions calculated on the lattice. With an eye towards lattice 
QCD attempts to describe nuclei, it is worth briefly noting the fundamental difference between 
lattice-measured correlation functions involving mesons and baryons. As an example, consider the 
noise to signal ratio of a correlation function involving n pion fields, where the small interaction is 
neglected. 



a{t) ^{A^-ADe-'^-^' , ^^^^ 



(0(0) ^/NAQ e-""'"' VN ' 

Here {d{t)) is the correlation function, a{t) is the variance and the A,- are amplitudes. It is note- 
worthy that in this ratio, the time dependence of the variance mirrors the time dependence of the 
correlator itself. One therefore concludes that correlators involving arbitrary numbers of pions have 
time-independent errors, as is indeed observed in lattice calculations. This of course renders the 
study of mesonic correlators quite pleasurable (from the statistical perspective). 

The baryons provide a more disturbing story; consider the noise to signal ratio for a proton 
correlation function: 



(0(0) VNAoe-'""' Vn 

Here the variance is dominated by the three-pion state rather than by the proton, and therefore the 
noise to signal ratio of the proton correlator grows exponentially with time. More generally, for a 
system of A nucleons, the noise to signal ratio behaves as 

(0(0) Vn ■ 

Therefore the situation worsens as one adds nucleons. Indeed one can show that the theoretical 



expectation from eq. |23| compares favorably with A = 2 data calculated by the NPLQCD collabo- 
ration [11]. 

The estimates given above, which follow from very general field theoretic arguments, naively 
indicate that nucleon and nuclear physics require exponentially more resources than meson physics 
to achieve the same level of accuracy. While these general results are somewhat discouraging to 
efforts to extract baryon interactions from the lattice, one should keep in mind that the arguments 
are rigorous only in the asymptotic, infinite time, limit. And, indeed, as we shall see below, there 
can be a region of intermediate times where these arguments do not apply, rendering the problem 
much less severe than these arguments suggest. 

3. Meson-meson interactions 

Here I will briefly review a calculation of the 1 = 2 pion-pion (tttt) scattering length, as it serves as 
a benchmark calculation with an accuracy that can only be aspired to at present for other systems. 
Due to the chiral symmetry of QCD, pion-pion nn scattering at low energies is the simplest and 
best-understood hadron-hadron scattering process. The scattering lengths for nn scattering in the 



s-wave are uniquely predicted at leading order in chiral perturbation theory {X-VT) \ 12] 



mna'^n = 0.1588 ; m^a'-j^ = -0.04537 , (3.1) 
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at the charged pion mass. While experiments do not provide stringent constraints on the scattering 
lengths, a determination of s-wave nn scattering lengths using the Roy equations has reached a 



remarkable level of precision [ ]13| , |14| ] : 



mno'-^ = 0.220 ± 0.005 ; ninJ-^ = -0.0444 ±0.0010 . (3.2) 



The Roy equations [ |15| ] use dispersion theory to relate scattering data at high energies to the scat- 
tering amplitude near threshold. At present, lattice QCD can compute nn scattering only in the 
1 = 2 channel as the 7 = channel contains disconnected diagrams. It is of course of great in- 
terest to compare the precise Roy equation predictions with lattice QCD calculations. Figure ^ 



summarizes theoretical and experimental constraints on the s-wave nn scattering lengths 014[]. It is 
clearly a strong-interaction process where theory has outpaced the very-challenging experimental 
measurements. 

The only existing fully-dynamical lattice QCD prediction of the I = 2 nn scattering length 
involves a mixed-action lattice QCD scheme of domain-wall valence quarks on a rooted staggered 
sea. Details of the lattice calculation can be found in Ref. ||^. The scattering length was computed 
at pion masses, niji ~ 290 MeV, 350 MeV, 490 MeV and 590 MeV, and at a single lattice spacing, 
b ~ 0.125 fm and lattice size L ^ 2.5 fm The physical value of the scattering length was 
obtained using two-flavor mixed-action X-?T which includes the effect of finite lattice-spacing 



artifacts to ^{m^b^) and 0'{b'*) [|16|]. The final result is: 



mna^nn = -0.04330 ± 0.00042 , (3.3) 




Figure 1: The state of threshold s-wave nn scattering. Noteworthy are the red ellipse from the Roy equation 
analysis and the grey band from the direct lattice QCD calculation of the 7 = 2 scattering length, as discussed 
in the text. 
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where the statistical and systematic uncertainties have been combined in quadrature. Notice that 
1% precision is claimed in this result. The agreement between this result and the Roy equation de- 
termination is a striking confirmation of the lattice methodology, and a powerful demonstration of 
the constraining power of chiral symmetry in the meson sector. It would clearly be of great interest 
to see other (fully-dynamical) lattice QCD calculations of the s-wave kk scattering lengths using 
different types of fermions. The K^K^ and n^K'^ scattering lengths have also been computed by 
the NPLQCD collaboration. I refer the interested reader to Ref. [IT]] for details. 



4. Meson-Baryon Interactions 

Pion-nucleon scattering has long been considered a paradigmatic process for the comparison of %- 
PT and experiment. To this day, controversy surrounds determinations of the pion-nucleon coupling 
constant and the pion-nucleon sigma term. The K n interaction is important for the description of 



kaon condensation in the interior of neutron stars [17], and meson-baryon interactions are essential 
input in determining the final-state interactions of various decays that are interesting for standard- 
model phenomenology (See Ref. [[l^] for an example). In determining baryon excited states on the 
lattice, it is clear that the energy levels that represent meson-baryon scattering on the finite- volume 
lattice must be resolved before progress can be made regarding the extraction of single-particle 
excitations. 

While pion-nucleon scattering is the best-studied meson-baryon process, both theoretically 
and experimentally, its determination on the lattice is computationally prohibitive since it involves 
annihilation diagrams. At present only a few limiting cases that involve these diagrams are being 
investigated [|I9|]. Combining the lowest-lying SU (3) meson and baryon octets, one can form five 
meson-baryon elastic scattering processes that do not involve annihilation diagrams: tt+S", 
K^p, K^n, and ^"S". (Note that ^"r+ has the same quantum numbers as tt+E".) Three of these 
processes involve kaons and therefore are, in principle, amenable to an SU (3) heavy-baryon ;t-PT 



(HB;^-PT) analysis [20] for extrapolation. The remaining two processes involve pions interacting 
with hyperons and therefore can be analyzed in conjunction with the kaon processes in SU (3) 
HB;t;-PT, or independently using SU (2) HBx-PT. Below we will review recent work in Ref. [§ 
which has carried out this analysis. 

The calculations of Ref. [||] were performed predominantly with the coarse MILC lattices 
with a lattice spacing of Z> ~ 0.125 fm, and a spatial extent of L ~ 2.5 fm. On these configurations, 
the strange quark was held fixed near its physical value while the degenerate light quarks were 
varied over a range of masses. (These are the same resources that were used to calculate the 71+ 71+ 
scattering length described above.) 

The scattering lengths of the five meson-baryon processes without annihilation have been cal- 
culated to &{m\,^) in SU{3) YiRx-^T and are given by 



1 



An TUn + mj: 



Iniji 2m 



m 



3 1 



J 71 



+ -^Ci + ) + 8/J123 (At ) -TJ 

Jk Jk 



(4.1) 



1 



m-^ 



An mjt + ms 



J 71 J 71 J 71 



. m-i 



(4.2) 
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Quantity LO (fm) NLO fit (fm) NLO (NNLO fit) (fm) NNLO (fm) 





-0.2294 


-0.208(01)(03) 


-0.117(06)(08) 


-0.197(06)(08) 


anS. 


-0.1158 


-0.105(01)(04) 


0.004(05)(11) 


-0.096(05)(12) 


aKp 


-0.3971 


-0.311(18)(44) 


0.292(35)(48) 


-0. 154(5 1)(63) 




-0.1986 


-0.143(10)(27) 


0.531(28)(68) 


0.128(42)(87) 




-0.4406 


-0.331(12)(31) 


0.324(39)(54) 


-0.127(57)(70) 



Table 1: SU{3) extrapolated scattering lengths. The first uncertainty in parentheses is statistical, and the 
second is the statistical and systematic uncertainty added in quadrature. 
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m 



■IK JK 



(4.4) 



1 



471 m/f + m-~, 



+ -#Ci + %,so(At) + 8/^123 (m)^ 



fi 



JK JK 



(4.5) 



where the C's and the /j's are low-energy constants (LECs) and the ^?^'s are loop functions given in 
Ref. [||]. It is clear that this is an overconstrained system. Many fitting strategies are possible, as 
discussed in Ref. [||]. We found it convenient to rewrite the chiral perturbation theory formulas as 
polynomial expansions. For instance, in the case of n^'ZP and K^n, we formed the objects: 



NLO 



NNLO 























m0 





= 1 - Coiwi0 



1 -Coim0 -8/ji(A;( 



(4.6) 



(4.7) 



for each of the processes. Notice that the left-hand sides of these equations are given entirely in 
terms of lattice-determined quantities, while the right-hand side provides a convenient polynomial 
fitting function. Figure ^ plots these functions for several of the processes. The shift of the value 
of r from NLO to NNLO is dependent on the renormalization scale jJ., and therefore with the 
choice jj. = A), one would expect this shift to be perturbative if the expansion is converging. The 
large shifts in F from NLO to NNLO are indicative of large loop corrections. The LECs fit to 
the lattice data are tabulated in Ref. [||]. While the NNLO LECs hi and hi23 appear to be of 
natural size, the NLO LECs Co and Coi are unnaturally large. The extrapolated values of the five 
scattering lengths are given in Table [I| While the TT+Z^ and n^E^ scattering lengths appear to be 
perturbative, the extrapolated kaon-baryon scattering lengths at NNLO deviate by at least 100% 
from the LO values. The seemingly inescapable conclusion is that the kaon-baryon scattering 
lengths are unstable against chiral corrections in the three-flavor chiral expansion, over the range 
of light-quark masses that we consider. 

Given the poor convergence seen in the three-flavor chiral expansion due to the large loop 
corrections, it is natural to consider the two-Havor theory with the strange quark integrated out. In 
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this way, TiL and n'E may be analyzed in an expansion in m^r. To ^{m\) in the two-flavor chiral 



expansion, one has [ ]23| ] 



1 



Ik niji + mz 



m 



-J2 + ^1^11+1.+ + -Ji^j^+Y,^ 



m-1 



K^z^ = j2^l + K^z^; (4.8) 



1 



471 rriji + 



ml ml 



K.E0 = y2^4 + K^Z0 , (4.9) 



where £4 is the LEC which governs the pion mass dependence of /^^ 1 13]. Note that the chiral logs 
have canceled, and in this form, valid to order in the chiral expansion, the scattering lengths 
have a simple polynomial dependence on m^j |23]. Figure || shows the 68% and 95% confidence 
interval error eUipses in the h-C plane for both 7r+Z+ and TT+S". Exploring the full 95% confidence 
interval error ellipse in the h-C plane yields 



= -0.197 ±0.017 fm ; 
a^+^o = -0.098 ±0.017 fm . 



(4.10) 
(4.11) 



These are the numbers that we quote as our best determinations of the pion-hyperon scattering 
lengths. In Figure ^ we plot the scattering length versus the pion mass. (Note that the bar denotes 
the scattering length rescaled by a kinematical factor 




71 Z SU(3) 



300 400 SOO 

(MeV) 

K*n SU(3) 




600 700 




Figure 2: Plots of Tnlo and Tnnlo versus the Goldstone masses for four meson-baryon processes in the 
SU(3) case. The line at F = 1 is the leading order curve, and dotted line is the physical meson mass. 
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012340123 
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Figure 3: The 68% (light) and 95% (dark) confidence interval error ellipses for fits for the 7r+E+ (left), 
and tt+E" (right) processes. 

a^.^. SU(2) a^.^o SU(2) 




700 
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Figure 4: a plots for the ;r+E+, and ;r+S processes versus the pion mass. The diagonal line is the leading 
order curve, and the dotted line is the physical pion mass. The innermost error bar is the statistical uncertainty 
and the outermost error bar is the statistical and systematic uncertainty added in quadrature. The filled bands 
are the fits to the LECs in the SU(2) case at NNLO. 



5. Baryon-Baryon Interactions 

The NPLQCD collaboration has performed the first full QCD calculations of nucleon-nucleon in- 
teractions [Q] and hyperon-nucleon |Q] interactions at low-energies with large pion masses. The 
nucleon-nucleon scattering lengths were found to be of natural size for the unphysically large val- 
ues of the quark masses at which these calculations were performed. In contrast, the physical 
values of the scattering lengths are unnaturally large compared to the range of the nucleon-nucleon 
interaction. This suggests a fine-tuning of the underlying parameters of QCD (the quark masses 
and QCD scale) and lattice calculations with quark masses much closer to the physical values are 
needed to reproduce the experimental values. By contrast, very little is known about the interac- 
tions between nucleons and hyperons from experiment, and lattice QCD calculations can provide 
the best determinations of the corresponding scattering parameters and hence determine the role of 
hyperons in neutron stars. 

The NPLQCD collaboration has undertaken an extensive exploration of the impact of high- 
statistics on one-, two- and three-baryon correlation functions on one ensemble of anisotropic 
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clover gauge configurations that are presently being generated by the Hadron Spectrum Collab- 
oration [^,10]. A total of ~300,000 sets of measurements were made using 1200 gauge configu- 
rations of size 20^ X 128 with an anisotropy parameter t, = bg/bt = 3.5, a spatial lattice spacing of 
bs = 0.1227 lb 0.0008 fm, and pion mass of Mj^ ~ 390 MeV. The ground state baryon masses (in 
lattice units) were extracted with uncertainties that are at or below the ~ 0.2%-level. An example 
of the precision that has been obtained for single baryon masses with this calculation is shown in 
Figure ||, where the generalized effective mass plot of the E baryon is shown. 
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Figure 5: Effective mass plot for the E baryon from a high-statistics calculation using anisotropic clover 
gauge configurations. 

With such statistics we were able to systematically explore the signal-to-noise issue in baryon 
systems, and found that, while signal/noise does degrades exponentially at large times, there is an 
intermediate time-interval for which it is time-independent, producing a noise-to-signal ratio per 
baryon that is essentially independent of baryon number. This is seen clearly in Figure ^ This 
feature greatly increases the number of nucleons that can be explored with lattice QCD for fixed 
computational resources. This was a significant finding and opens up the possibility of calculations 
with four and more baryons. In work that is currently in progress, we have obtained a number of 
phase-shifts in the two-baryon sector with > 5 sigma-level significance. Figure |7| gives the energy 
shifts of various processes with associated statistical and systematic errors added in quadrature. 



6. Conclusion 

Lattice QCD calculations of two- and three-body interactions of pions and kaons are now a pre- 
cision science (for those channels that do not involve disconnected diagrams). I have reviewed 
the first fully-dynamical lattice QCD calculation of meson-baryon scattering. An analysis of the 
scattering lengths of these two-body systems using HB;^PT has led to the conclusion that the three- 
flavor chiral expansion does not converge over the range of light quark masses that are investigated. 
By contrast, the n^Yf^ and 71+ S" scattering lengths appear to have a well-controlled chiral expan- 
sion in two-flavor Y{&xVY. Our results, an+Y,+ = -0.197 ± 0.017 fm, and a^+^o = -0.098 ± 0.017 
fm, deviate from the LO (current algebra) predictions at the one- and two-sigma level, respectively. 
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Clearly a milestone for this area of research would be to see a definitive signal for nuclear 
physics. In the case of the NN interaction, in addition to signal/noise issues, one must face a fine- 
tuned system that requires a non-perturbative effective field theory description. While this poses a 
tremendous challenge, this year has seen a great deal of progress in the NPLQCD collaboration's 
state goal to compute nuclear properties using lattice QCD. During the last year we performed 
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Figure 6: Effective mass plots for various one-, two- and three-body correlation functions calculated using 
anisotropic clover gauge configurations. There is an intermediate region of time slices for which signal/noise 
is time independent. 
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Figure 7: Energy splittings for various two- and three-body processes calculated using anisotropic clover 
gauge configurations. 
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a high-statistics calculation of one-, two and for the first time — three-baryon systems 

which represent a "jump" by an order of magnitude in the data volume produced in such calcu- 
lations. The remarkable conclusion that the signal/noise problem is not as severe as previously 
thought promises to hasten a golden age of exploration for nuclear physics using lattice QCD. 
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